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Cauchy problem for integrable discrete 
equations on quad-graphs 

Dedicated to S.P. Novikov 
on his 65 birthday 


V.E. Adler* A.P. Veselov^ 


Abstract 

Initial value problems for the integrable discrete equations on quad- 
graphs are investigated. A geometric criterion of the correctness of 
such a problem is found. The effects of the interaction of the solutions 
with the localized defects in the regular square lattice are discussed for 
the discrete potential KdV and linear wave equations. The examples 
of kinks and solitons on various quad-graphs, including quasiperiodic 
tilings, are presented. 


1 Introduction 

The discrete potential KdV (dKdV) equation 

('^m-|-l,n+l ^m,n+l) Olm fdni TTl^Tl G Z (1) 

can be considered as a simplest representative of the integrable nonlinear 
discrete equations in two dimensions. This held is widely studied and many 
other examples can be found in the literature, see e.g. [1-7]. However, we 
restrict ourselves to this simplest model, since our aim here is to analyze the 
generalizations of the other kind. 

*Landau Institute for Theoretical Physics, 12 Institutsky pr., 142432 Chernogolovka, 
Russia. E-Mail: adler@itp.ac.ru 

^Loughborough University, Loughborough, Leicestershire LEll 3TU, UK and Landau 
Institute for Theoretical Physics, Russia. E-Mail: A.P.Veselov01boro.ac.uk 


1 



Namely, we will discuss the role of the support of the discrete equations, 
that is, roughly speaking, the set where the independent variables live. Of 
course, the most natural way to discretize 2-dimensional PDE is to consider 
equations on the regular square grid 7? as above, but recent results [8 - 15] 
demonstrate that possibly more general synonym of ‘2D’ in the discrete case 
is ‘planar graphs’. 

It should be noted that non-standard planar graphs were considered in 
Mathematical Physics already quite a while ago (see e.g. Korepin’s works 
nisi, [T^ where the solvable spin models on the quasi-crystallic tilings were 
investigated). A systematic theory of the linear difference operators on the 
graphs in relation to soliton theory has been initiated by S.P. Novikov (see 
|13 - 17]). 

We start in Section with the necessary information about integrable 
equations on quad-graphs, which are planar graphs with quadrilateral faces. 
In this case the integrability can be understood as the so called 3D consistency 
property, introduced and studied in [9 - 12]. The discrete potential KdV 
equation on quad-graph reads locally exactly as on the square lattice: 


(ni2 - n)(ni - ua) = «! - 02, (2) 

but now it is assumed that the helds v are assigned to the vertices of the 
quad-graph and the parameters a are assigned to the edges, as shown on the 
Fig.0 below, and the relation (^ must be fulhlled for each face of the graph. 

Section ^ is devoted to the general discussion of the possible settings of 
Cauchy, or initial value problem (IVP) on quad-graphs. At the hrst sight 
the question about the correctness of IVP is not related to integrability but 
in fact it does very much (see e.g. Fig. |T^, below). The 3D consistency 
property plays an essential role here. 

Section ^ contains the main result of the paper: a criterion for the ex¬ 
istence and uniqueness of the solution of Cauchy problem for an integrable 
equation on a quad-graph (Theorem 

In section 5 we discuss the effect of interaction of the solutions of the 
dKdV equation with the localized defects in the regular square lattice. We 
introduce a notion of the weak defect and show that such a defect does not 
affect the dynamics outside it. 

Several examples of the kink/solitons solutions on various quad-graphs, 
including quasiperiodic tilings are presented in the section ^ 

In the last section we discuss the solutions of the linear discrete wave 
equation 

1^12 - '^1 - 1^2 + W = 0 

on quad-graphs. Surprisingly enough, in contrast to the non-linear case these 
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solutions have a non-trivial interaction with weak defects, although the effect 
is limited (see Theorem ^). 


2 Quad-graph equations 


We start with the results of the paper [|I^ where all the integrable cases 


were classified among the certain class of discrete equations on quad-graphs. 
The basic building block of such equation is an equation on the elementary 
quadrilateral of the form 


(5(n,ni,n2,t’i2;ai,a2) = 0, (3) 

where the held variables n G C are assigned to the vertices and the parameters 
a G C are assigned to the edges of quadrilateral, as shown on the Fig.|l|. The 
following properties are assumed to be satished: 

1) the parameters a on the opposite edges of any face are equal; 

2) the equation (^) is invariant under the group of the square symme¬ 
tries: 

Q{v, vi, V 2 , V 12 ; «!, 0 : 2 ) = £Q{v, V 2 , Vi, Vu] 02 , ai) 

= aQ(vi,v,vi 2 ,V 2 ;ai,a 2 ), e,a = ±l; (4) 



ai 

Figure 1: An elementary quadrilateral and its symmetry 


3) the function Q is a polynomial of degree one with respect to each held 
variable. 

The hrst property means that the value of parameter is constant along 
any strip (or, characteristic), which is a sequence of quadrilaterals adjacent 
by the opposite edges. We will see in the next sections that this notion is 
very important in the theory of quad-graph equations and plays the role 
similar to the characteristics for the second-order hyperbolic eqautions. In 
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the particular case of dKdV equation on the square lattice (|^ we have two 
families of the parameters and jSn assigned to the vertical and horizontal 
straight strips, but in the general case the strips may bend and intersect in 
a very complicated way (for instance, self-crossing, self-tangent and closed 
strips are possible). Some characteristics are shown at the pictures below by 
dashed lines crossing the adjacent edges of the faces in the strip. 

The second property allows us to define equation on each face indepen¬ 
dently on its position in the quad-graph. 

The third property means that the equation (|^) can be solved with respect 
to a field at any vertex as linear fractional expression of other three fields. 
This expression could be singular for some special values of the helds, which 
is not essential for many types of quad-graphs, since we are interested only 
in the generic initial data. However, this becomes important for the graphs 
containing self-crossing strips, as we will see in Example 5. 

The notion of integrability laid in the basis of the classification in 
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is the 3-dimensional consistency |P, which means that the equation 

d^) may be consistently embedded into a 3-dimensional lattice, so that the 
similar equations hold for all six faces of any elementary cube, as on Fig.|^. 
It should be mentioned that this property is in a close relation to the set- 
theoretical solutions of the Yang-Baxter equation [^, or Yang-Baxter maps 


(see the discussion of this in the Conclusions of 


V2, 



V23 


V2 


Figure 2: Three-dimensional consistency 


To describe more precisely what does the 3D consistency mean, consider 
the Cauchy problem with the initial data v,vi, 02 , 03 . The equations on the 
rear faces 

Q{o, Oi, Oj, Oij] ai, aj) = 0, I < i < j < 3 (5) 
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allow one to determine nniqnely the valnes ni 2 ,'^ 23 - After that one has 
three different eqnations for ^ 123 , coming from the front faces, and consistency 
means that all three valnes thns obtained for ^123 coincide. For example, in 
the discrete KdV case we have the formnla 

^ _ (ftl - a2)^’l^^2 + (Q3 - ai)V'iVl + (0^2 - Q3)^2t^3 

(q! 3 - tt2)vi + (ttl - a3)^'2 + (tt2 - ai)vz 


independently on the order of calcnlations. 

It should be noted that the variable v does not appear in this expression. 
This property is very important in many aspects, in particulary in the clas- 
sihcation problem which was solved in only under this assumption. The 


role of this property for the Cauchy problems is not very clear. Although we 
will not use it in our analysis of the IVP below we will see that some results 
for the linear wave equation (which satishes all the properties l)-3) but not 
the n-independence) and for the equations from the list will be different. 
Of course, linear equation is exceptional from many points of view, so it does 
not mean that this is related to this particular property. 

The following Theorem is the immediate consequence of 3D consistency. 
It defines the Backlund transformation for quad-graph equation, and we will 
make use of it in the next Section. 


T 



Figure 3: Backlund transformation 


Theorem 1. Let variables Vi, i G D(F) satisfy 3D consistent quad-graph 
equation & on the quad-graph F. Let us fix some vertex zq and choose 
arbitrarily the values Vi^ and a. Then the formulae 

Q{vi,Vj,Vi,Vj;aij,a) = 0, {i,j)eE{T) (7) 

define uniquely new solution Vi of equation (Qj. 
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Proof. Let us consider some path from io to the given vertex in- Along this 
path we have the sequence of equations of the form 

which dehne recurrently ..., The fact that the different choice of the 
path gives the same result and that the constructed Vi solve (|) follows from 
the 3D consistency. □ 

The discrete zero curvature representation for equation (^) 

= Lkj'Lji^i, (j, /c), {i,j ), (j , k) G T'(r) (8) 

where Lm,n = L{vn, Vm', /) is 2 X 2 matrix can be obtained by interpreta¬ 
tion of the new solution on the upper copy of T as the -^-function i.Q.0- 
To do this, the formula (0) at a = / (‘spectral parameter’) should be rewritten 
in the form of the discrete Riccati equation 


Vj = R{vi), {i,j)eE{T) 


where the function R is linear fractional due to the property 3), with the 
coefficients depending on Vi, vj, aij, I, and the latter equation is linearized by 
the substitution v = fj/cf). In the discrete KdV case this procedure brings to 
the matrices 


^j,i ~ 


ViVj -|- a 




( 9 ) 


3 Initial value problems 

Let us consider the IVP (or Cauchy) problem for a discrete equation on a 
quad-graph T with initial data on a connected subgraph P. In most of the 
examples we consider P is a simple path which is a connected sequence of 
the edges of the graph without self-intersections. 

Slightly abusing terminology we say that such a Cauchy problem is cor¬ 
rect if for the generic initial data the solution exists and unique (without 
investigation of the stability question). For hnite graphs this is certainly 
justihed since the stability is automatic then. 

Let us hrst compare the number of equations and unknowns for the given 
quad-graph T. We assume that T is simply connected and hnite, and contains 
F faces and V vertices, from which 14 ones belong to the boundary. Now 
apply the double counting of the angles in all faces (we assume for simplicity 
that each edge is a straight line, which, of course, is not necessary; the Euler 
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formula provides the universal method). Since the sum of angles in each 
quadrilateral is equal to 27r, hence the total sum is equal to 2nF. On the 
other hand, summation over the vertices gives 2n{V — 14) for the interior 
nodes plus 7r(14 — 2) for the boundary, so that the following equality holds: 

F = V -Wb-l. 

2 

This means that in order to balance the number of equations and unknowns 
one has to to assign initial data in some ^14 + 1 vertices. These vertices are 
not necessarily he on the boundary, as the simplest example of a staircase in 
rectangular domain of the square lattice demonstrates (see Fig.^; here and 
everywhere initial data are marked by solid line). 


Figure 4: Possible choices of initial data on the square grid 


Remark. The same calculation in the case of planar graph with /c-lateral 
faces yields the formula {k — 2)F = 2V — Vb — 2. It means, under the natural 
assumption that Vb «V for large domains, that graphs with triangular faces 
lead to the overdetermined systems, while systems on graphs with k > A are 
underdetermined. Of course, the balance between equations and unknowns 
can be restored in many ways. For example, one can consider systems with 
variable number of equations per face. In this way some integrable models 
were obtained on the triangular and hexagonal lattices, see e.g. |^. Other 
possibility is to consider graphs with various types of the faces. 


Of course, the above consideration is too rough and cannot answer how to 
choose the correct settings of Cauchy problem even in the simplest situations. 
In order to get some experience, let us consider several examples assuming 
that equation (|[) satishes the properties l)-3) from the previous Section, but 
is not necessarily integrable at hrst. 

A priori, several cases are possible: the IVP may be correct, or it may 
be underdetermined or overdetermined, and the calculation scheme for the 
correct IVP may be explicit (need only solving of one basic equation (^) 
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at each step) or implicit (need to solve the systems of such equations). It 
turns out that all these possibilities can be easily realized on very simple 
quad-graphs. 




Figure 5: Correct and overdetermined Cauchy problems 


Example 1. It is easy to see that Cauchy problem shown on the Fig. ^ 
is correct, that is all helds are uniquely dehned for arbitrary initial data 
(except for possible singularities which may arize in dKdV equation due to 
zero denominator). In contrast, the IVP on the right picture is not correct. 
Indeed, the values n, Ui, V2 are uniquely dehned, but for the values tc, V12 we 
have the system of two equations (use symmetry (|D): 

Q{y, ui, w, ni 2 ; a, j3) = 0, Q{y, V 2 , w, vu, a,/3) = 0 

which gives the contradiction in the generic case vi ^ V 2 - 

Example 2. Both IVP shown on the Fig.^ are correct. For the left one 
this is obvious; for the right one we have the system of three equations on 
the variables m, w,vi 2 : 

Q{v, ui, M, vu] a, /3) = 0, Q{v, w, m, vu, a, /3) = 0, Q{v, w, V 2 , vu, a, /3) = 0. 

This system implies u = V 2 , w = Vi and for V 12 one has the same equation 
Q{v,vi,V 2 ,vi 2 ', C(, P) = 0 as in the case of lattice without defect. In other 
words, the defect shown on these pictures does not affect the dynamics. We 
shall call defects with such property transparent. 

Analogously, the graph shown on the Fig. ^ gives explicit and implicit 
schemes depending on the choice of initial data. 


































Figure 6: Correct Cauchy problems with explicit and implicit calculation 
scheme on transparent defect 




Figure 7: Correct IVP with explicit (a) and implicit (b) schemes 


Example 3. For the quad-graph shown on the Fig. ^ the hrst two IVP are 
incorrect: for a) case one hnds successfully the values v inside the square 
bounded by initial data and the dashed characteristic, but encounters cotra- 
diction outside it; for b) case the same square remains undetermined. 

Example 4. Several examples of correct IVP on different quad-graphs are 
shown on the Fig.^. Inspired by these examples and analogies with continu¬ 
ous case one can make a conjecture that an IVP with initial data given on a 
simple path P is correct iff each characteristic (or strip) in P meets P exactly 
at one edge. 

We hnish this Section with two examples which show that for the general 
equation this conjecture is far from being true and emphasize the role of the 
3D consistency condition in this problem. 
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Figure 8 : Overdetermined, underdetermined and correct Cauchy problems 


Example 5. On the quad-graph on the Fig.|^ the vertical characteristics 
come back and cross the initial data again, while the horizontal characteristics 
do not cross the initial data at all. Nevertheless, this IVP is correct in 
general case, by means of implicit calculation scheme which can be described 
as follows: 

1 ) the implicit step: unknowns U 5 , Ue, Vj, Vg are found as functions on Ui, V 2 , 
■^^ 3 ) 'V 4 by solving the system 


Q{vi,V 2 , U 5 , Vq] a, 7 ) = 0 , Q{v 2 , V 3 , ve, vr, /3, 7 ) = 0 , 

Q(v3,t'4,t'7,t'8;a,7) = 0, Q{v&,v7,v^,vg-,j3,a) = Q. (10) 


2 ) next, all other variables on the bottom strip are computed explicitly; 

3) now we are in the same situation and can repeat procedure. 


Of course, the question is if the above system is degenerate or not. There 
is no reason for it to be degenerate for generic Q, but it turns out to be 
degenerate for 3D consistent equations. Indeed, in this case we can con¬ 
sider four quadrilaterals corresponding to the system (p!0[) , as the faces of 3D 
cube, and therefore the initial data Ui, U 2 ,'^ 3 , 1^4 must be related by equation 
Q{vi,V 2 ,Vi,vg]a, P) = 0 in order to provide compatibility of this system. 
Moreover, even if this constraint is fulhlled, the values V 5 ,vq,vt,vs cannot 
be found uniquely. So, in the integrable case this IVP can be considered as 
overdetermined and underdetermined simultaneously, in accordance to our 
conjecture. 


Example 6. The self-crossing strip on the Fig. 0 shows that we need to 
make some stipulation in the integrable case. In this example the path with 
initial data U 2 ,U 3 ,U 4 ,U 5 crosses 3 characteristics; one computes consequently 
Vi,Vq,Vj,V 8 and Vg = automatically due to the 3D consistency. At this 
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Figure 9: Examples of correct Cauchy problems 


point the dKdV turns into (uio—'y 4 )('i^ 9 —'^ 5 ) = 0, so that the value Uio remains 
arbitrary. However, in this case we can eliminate the self-intersecting strips 
from consideration by splitting them into the pairs of independent strips (one 
of them closed) by identification of the vertices in which the fields coincide, as 
shown on the Fig. 0 It should be noted that such splitting of the equation 
on the quadrilateral with the pair of equal parameters: 


Q{v,vi,V2,vu;a,a) = 0 {vu - v){v2 - vi) = 0 


takes place for all equations from the list obtained in the paper 


4 Existence and uniqueness theorem 

Now we are ready to formulate our main result. 

Theorem 2. Let equation (Q) he 3D consistent, and let P he a simple path in 
a finite simply connected planar quad-graph F without self-inters eating strips. 
Consider Cauchy problem for this equation with generic initial data on the 
path P. 

1) If each characteristic in F intersects P exactly at one edge then the IVP 
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Figure 10: Correct/overdetermined Cauchy problem (general/integrable 
case). 


10 I 9 I 8 
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Figure 11: Nonuniqueness in uio 


is correct, that is the solution exists and unique. 

2) If some characteristic intersects P more than once then the IVP is overde¬ 
termined (no solution exists for generic initial data). 

3) If some characteristic does not intersect P then the IVP is underdeter¬ 
mined (if a solution exists, it is not unique). 

We should clarify that when we consider the intersections of characteris¬ 
tics with a path we count only the internal edges of the characteristics. For 
example the intersection of a vertical characteristic with a staircase is one 
edge (but not two). 

We will prove this theorem by use of some transformations of quad-graphs, 
based on the notion of Backlund transformation (see Theorem |I|). However, 
the existence of solution can be proved also by use of quite different idea of 
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Figure 12: Elimination of self-intersecting strip 


immersion of F into the multidimensional unit cube. 

Proof of existence. Let the conditions of the part 1) of the Theorem are 
fulhlled. Let us construct an immersion of F into the iV-dimensional unit 
cube, where N is the length of P. To do this, assign the shift operators 
Tj, i = 1,... ,N to the edges of P. By dehnition, Tj acts on iV-dimensional 
vectors by formula Tfx) = x + Ci (mod 2) where e* is the vector with 1 at 
Tth place and 0 at the others. 

After this, let us distribute these operators at all edges of F, accordingly to 
the same rule as for the parameters a, i.e. one operator per strip. Obviously, 
this is possible in unique way due to the condition 1). Now, let us choose 
the first vertex uq of P as the origin with coordinates (0 ,..., 0) and dehne 
the coordinates of any other vertex v as the result of applying shift operators 
along any path from vo to v. The result does not depend on the choice 
of P^. Indeed, any closed path has an even intersection index with any 
characteristic and therefore is equivalent to identity operator. 

So, we have found an immersion of F into iV-dimensional unit cube, and 
each face of F now is just a 2-dimensional face in this cube. Initial values are 
given on the path 

P=((0,0,...,0), (1,0,..., 0), (1,1,...,0), ..., (1,1,...,!)) 

and allow to calculate the values in each vertex of the cube. (It is not difficult 
to explain how to obtain the values for the coordinate vectors e^, and so on.) 
Of course, 3D consistency is necessary here in order to justify the correctness 
of the procedure. In this way, we have found some solution of the Cauchy 
problem. □ 

Remark. It is worthy to note that all the helds at the vertices of the multi¬ 
dimensional cube can be calculated recursively using only one equation each 
time, so that on the cube the integration scheme becomes explicit. However, 
this algorithm of solving our Cauchy problem is not very effective, since the 
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number 2^ of the vertices in the A^-dimensional cube grows exponentially in 
N, while the number of vertices in F (and therefore the size of the algebraic 
system we are solving) is about We would like to mention also that a 
general question of embedding of a quad-graph into cubic lattice was investi¬ 
gated in details by Dolbilin et al |23]. The absence of self-intersecting strips 


is crucial for the existence of such embedding. 


The general proof presented below is based on the simple transformations 
of erasing and inserting of a strip. These transformations bring, roughly 
speaking, to equivalent IVP and can be described as follows. 



Figure 13: Erasing and inserting of a strip 


Let us consider some strip C (with parameter a). It divides T into two 
subgraphs, Ti and r 2 . Let us denote Pi the corresponding borders of the 
strip and let ii G Pi and *2 G P 2 be a pair of neighbour vertices separated 
by C. Now, let us apply BT (0) to the part r 2 of the quad-graph, taking *2 
as the seed point, and Vi^ = This will give us some new solution in r 2 
such that the variables v along the border P 2 coincide with the corresponding 
variables v along Pi. Therefore we can construct the common solution on the 
graph T obtained from T by removing of the strip C. This solution coincides 
with the old one on Ti. Of course, some information is lost after such an 
operation and in order to make it invertible one have to remember the value 
Vi2- 

On the Fig.|^ these operations are shown for the case of simple strip, 
however the closed or self-tangent strips are allowed as well. For the self¬ 
tangent strips the 3D consistency property implies the coincidence of the 
helds in some vertices, and these vertices can be merged during the above 
procedure (Fig. p!^). 
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Figure 14: Case of self-tangent strip 


Proof of necessity. 2) Let some characteristic C intersects the path P 
with initial data twice. Let us enumerate the vertices of P as shown on the 
Fig. |T^ so that the edges of intersection are (1,2) and {n — l,n). 



Fi 




Fi 


Figure 15: Constraint on the initial data 


Suppose that IVP is correct, then solution of equation exists on 
the subgraph Fi bounded by C and P. Let us look at this picture as 3- 
dimensional one, assuming that the exterior border of the strip C lies at the 
next floor, in particular the vertex 1 is over 2 and n is over n — 1. Accord¬ 
ingly to the Theorem |l| the value can be dehned by solving the sequence 
of equation (|^ along the path (2 ,..., n — 1), that is is some function on 
Ui,..., Vn-i- This means that the initial data are not free, that is the IVP is 
overdetermined. 

2) Suppose that some strip C never intersect the path P with initial data. 
Let US cut off this strip, as shown on the Fig. |T^ (we assume for dehneteness 
that P lies in the subgraph Fi). If a solution of the IVP exists in P then 
this procedure gives ns a solution with the same initial data on the reduced 
graph P. 

Next, let US apply the inverse procedure to this solution, that is consider 
BT on P 2 , taking some arbitrary value for Vi^. This gives ns some solution 
on the graph P which coincides with the original solution on Pi but not on 
P 2 . This means that the initial data are not sufficient to provide uniqueness 
of the solution. □ 
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Proof of uniqueness. This can be obtained by seqnential erasing of the 
strips, as explained above, nntil the qnad-graph is exhansted. □ 

Remark. Althongh we assnmed that P is a simple path (which is probably 
the most natnral case), it follows from the proof that the resnlt is trne for 
any connected snbgraph. 

We will call a connected snbgraph P Cauchy subgraph if it has the prop¬ 
erty that each characteristic in T meets P exactly at one edge. 

A natnral qnestion is whether snch a snbgraph does always exist or not. 
The answer in general is negative. The following pictnre gives of an example 
of a regnlar lattice on which no correct IVP is possible, at least with initial 
data on the connected snbgraphs. It is obtained from the kagome lattice 
by the standard process of merging with the dnal graph (which is the qnad- 
graph shown on the Fig.||d). For the eqnations on qnad-graphs of snch type 
the qnestion arise how to describe the set of general solntions. 



Fignre 16: Example of transitive lattice withont connected Canchy snb¬ 
graphs: all characteristics are closed 


5 Interaction of the solutions with the local¬ 
ized defects 

Consider the homogeneons dKdV eqnation 

i^^m+l,n+l '^m,n)(Cm+l,n '^m,n+l') ^ /3, m, H ^ 'Z (H) 

on the regnlar sqnare lattice, assnming that all horizontal strips carry one 
and the same parameter a and analogonsly for the vertical strips. 
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Let us replace now some m x n rectangle by a finite quad-graph D (with 
the same vertices on the boundary). What we will get is a regular square 
lattice with localized defect. Many examples of such localized defects can be 
found in the previous section, see e.g. Figures H) B 0) Outside of the 
defect we have a regular picture of horizontal and vertical strips. 

We will say that the defect is weak if all the characteristics entering the 
defect leave it in the same direction but possibly in a different order. One 
can see the examples of weak defects on the Figures On the Fig. |] 

we have a non-weak defect since a pair of characteristics change their types. 

It is clear that the homogeneous dKdV equation can be naturally defined 
on a regular lattice with a weak defect. Consider for the corresponding 
equation a Cauchy problem with the initial data outside the defect (e.g. on 
the staircases) just like for the usual square lattice. By the results of the 
previous Section for generic initial data it has a unique solution. Let us 
compare this solution with the solution of the same Cauchy problem but on 
the regular square lattice without defects. 

The result turns out to be surprising: outside of the defect they will be 
the same. 

Theorem 3. For the homogeneous dKdV equation a weak defect does not 
affect the generic solution outside the defect. 

In other words for homogeneous dKdV solutions the weak defects are 
transparent. This is especially remarkable since for the discrete wave equa¬ 
tion this is not true (see below Section 7). 

The proof is based on the existence of discrete zero-curvature representa¬ 
tion and follows from the following general considerations. 

Let P be a closed path in a simply connected quad-graph F then for 
any solution of dKdV equation on F the ordered product of the transition 
matrices (|^) is equal to a scalar matrix. More precisely the formula 

p c 

is valid where the product in the right hand side is taken over all character¬ 
istics in F, ttc is the parameter corresponding to the strip C and the index 
ind(C', P) is the total number of intersections of C and P (without signs). 
This formula can be easily proved from the property (H) and the relation 

k-'ijLji (/ 

Equivalently, if Pi and P 2 are two paths from zq to p then 

Pi c P2 
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Theorem 4. Let the paths Pi and P 2 have the common starting and end 
points and for any characteristic C in T we have ind(C, Pi) = ind(C, P 2 ) < 1. 
Then for any generic solution of dKdV equation on V the fields v on P 2 can 
be uniquely recovered from the values on Pi. 

Proof. Let us consider the matrix product L{1) = Lji = (“ ^). One can 
prove that for a generic solution rankL(ac') = 1 for each strip C intersecting 
Pi (for some special solutions this can fail: for example if we have 0^*2 = 0^2213 
and Ujg = Up for three consequent vertices in the path, then . = 

0 ). 

Now the result follows from the uniqueness of the refactorization of the 
same matrix as L{1) = Lji along the P 2 = ■ ■ ■ )• It is enough to 

show that there exists a unique matrix of the form 


Lilio 


-1 Vi 


such that L{1) = LLi^^i^ and the matrix L is polynomial in 1. Indeed, Vi^ is 
uniquely defined from the condition kerL(Q;jpp) ~ Since the polynomi¬ 
als aup -|- h and cvi^ + d are divisible by I — ajgjj the matrix 

1 f-avi^ - b (cTjop - l)a Vi^^aVi^ + bf 
. -cVi, - d 


L = 


I 


Oti. 


i ^7.1 


- l)c P ViQ^cVi^ Pd)^ 

is polynomial. Continuing in the same way we will reconstruct the solution 
on the path P 2 . □ 

Now to derive the Theorem one should consider two sides of the rect¬ 
angle D directed towards the Cauchy subgraph as Pi and other two sides as 

P2- 

Notice that in general the Cauchy problem with the initial data on Pi 
may be underdetermined, since some characteristics may not intersect it. For 
example, some closed strips may occur inside the subgraph bounded by Pi 
and P 2 . Nevertheless, the use of transition matrices Lji will allow us to skip 
over such regions of underdeterminancy. 

For the correct IVP the Theorem ^ gives an integration scheme which 
allows to obtain the solution of dKdV equation by refactorization of the L 
matrices product along the path of initial data. 


6 Kinks and solitons on quad-graphs 

In this section we consider some explicit solutions of dKdV equation rather 
than solutions of Cauchy problem. The hint is to search for solution com¬ 
patible with the continuous dynamics. The dKdV equation on lattice 
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d^) originates from the nonlinear superposition principle for the Backlund 
transformation 

Vx + V^ = {v- vf - I (12) 

of the potential KdV equation Vt = v^xx — Usage of equation o allows 
to obtain the prohle of the kink solution, and then to hnd algebraically the 
dependence of the phase on m and n (the dependence on t can be easily 
recovered, but now this is unessential). 

Let us start from the the simplest situation of constant parameters am, = 
a, (5n = P- In this case the linear seed solution is obvious 

= o,m + bn + px + q, — h^ = a — {3. (13) 

Applying the Backlund transformation (0 one obtains 

Vm,n = dm + bn + px + q — kta.nh.{kx + (f>m,n), I = 2p — 

Let us search for the dependence of (j) on m,n in the form 

(t>m,n = pm+ nn + ^ 


and denote 

Then substitution into (P^d) yields 


(^ + b + k 


X-l 
X + 1 


MNX +1 ) 


— b + k 


NX - 1 
NX + 1 


MX + l) 


2 1,2 
= a —b 


Collecting together and splitting over X give the relations between the pa¬ 
rameters which can written in the form 

atanh/i = fotanhi/ = k, 

so that hnally 1 -kink solution is given by the formulae 

Vm,n = dm + bn + px + q — ktanh.{kx - 1 - pm + nn + ^), (14) 


where 


2 12 n 1 a + k 

a -b =a-f3, jj,=-log --, 

2 a — k 


n = 



b + k 
b-k' 


Obviously, one has to choose \k\ < max{|a|, | 6 |} in order to obtain real /r, n. 
The dependence on the parameter x can be neglected, but we may also make 
use of it in order to obtain the discrete dynamics of the soliton solution Um,n = 
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—2dx(vm,n) which corresponds to the usual KdV equation ut = Uxxx + GuUx 
(see Fig. [T7|) . 

One can check using the formula (|^ that the corresponding helds u 
satisfy the following discrete KdV equation 

Vui 2 + U2-2a - Vui + u-2a = uu + ui -2(3 - \/ U 2 + u - 2(3 (15) 

or, equivalently, 

-(ni2 - uf{ui - U2Y - {a - ( 3 ){ui 2 - u){ui - U2){u + ui + U2 + M12) 

+ (a - (3Y{{ui 2 - uY + (mi - U 2 Y) + 2 (q!^ - (3'^){ui2 - u){ui - U 2 ) = 0 . 

(16) 

In contrast to the discrete potential KdV equation this equation is quadratic 


with respect to each field, and thus is not covered by the classification |12 


u 


Figure 17: Soliton at a = 1, 6 = 2, fc = 0.3, p = q = x = G^^ = —7.5 

Two-kink solution can be obtained via NSP for Backlund transformation, 
which brings to the formula 

^2 _ 

.y(12) ^ .y(O)_ ^2 ^1 _ 

k 2 tanh(A; 2 a; -|- 02 ) — ki tanh(fcia: -|- 0i) 
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where the seed solution and phases 02 are dehned as above. In order 
to obtain prohle without singularity in x one may choose k 2 > ki > 0 and 
make a shift 02 02 + y- la this way the general multikink solution can 

be constructed. Notice the different effect of the Backlund transform which 
may result either in the dressing of the solution or just in the phase shifts 
m—+ n ^ n + 1. 

The seed solution (^) is easily generalized for the case of variable pa¬ 
rameters: 

^m,n T T T Qi (^771+1 ®m) ipn+1 ^n) 0n 

(Notice that this is more general ansatz even in the case of constant a,/?.) 
Proceeding as above one obtains the kink moving with the variable velocity: 

Vm,n = ttm + bn+px + q- ktunh^kx + Pm + J^n + 0) (l^) 

where 


(^m+i (^m) tanh(^m_|_i Pm) (bn+1 ^n) tanh(i^ji_|_i l^n) k^ 


that is 

1 , ®m+l 0,m T k 

pm+1 pm TT lo§ T 1 

^ O^m+l Ojm ^ 

On the Fig. 1^ the graph of the soliton u 
the values 


1 , bn+i — bn + k 

-^n = - log 7-7- 

2 ^n+l bn k 


= —2vx is plotted corresponding to 


ao = bo = po = uo = P = q = X = 0 , 
am - ttm-l = a/tU, bn - 6n-l = k = 0.3, ^ = -3. 


Moreover, the formula (^) can be easily generalized for the case of mul¬ 
tidimensional lattice in which any 2D sublattice is governed by dKdV 
equation: 


Vn^,...,n^ = al^ + --- + a‘l^+px + q- kia.n\i{kx + + ■ ■ ■ + 0, (18) 

where 



<Y- 


«n = Y 


'n+l 


= 


«+! - O tanh(z/;+i 

Ijp + 

2 



k, 
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Figure 18: Bended soliton 


This observation immediately gives us the kink solution for the quadgraphs 
which can be embedded into lattice. For instance, the graph shown on 
the Fig. Id can be obtained from the cubic lattice by taking all its faces 
which intersect the plane ni + n 2 + = 0, so that the kink solution for 
this graph is given by the formula (^) at d = 3 restricted on the values 
n-i + n 2 + ns G { — 1,0,1}. In the case of the linear phase = z/*nj the 
corresponding soliton propagates in the direction — z/^, — z/^) 

(see Fig.|T9|). 

Some subclass of the quasiperiodic tilings can be obtained if one choose 
the other planes and use de Bruijn’s projection method (see e.g. [^). The 
formula ([I8| ) works for such quad-graphs as well. On the Fig. ^ the same 
3D solution is presented, restricted on the plane ni -|- \/2n2 + = 0. 

Two more hgures Fig.^ and present the 3-dimensional solution re¬ 
stricted on the coordinate quadrants. 

An interesting problem is to investigate the interaction of the kink so¬ 
lutions with the localized defects in the square lattice. According to the 
Theorem a weak defect in the homogeneous lattice does not affect solu¬ 
tions, so that the ordinary kink (|^ does not change. In contrast, it can 
be shown that the bended kink (pTf) may acquire some phase shift after the 
interaction with such a defect. 
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Figure 19: Soliton on the regular quadgraph (oi = 0.5, 02 = 2.5, 03 = 3, 
k = 0.3, p = q = x = ^ = Qi) 


23 






Figure 20: Soliton on the quasiregular quadgraph 
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Figure 21: Soliton on three quadrants n* > 0 (oi = 1, 02 = 1.5, 03 = 2, 
k = 0.3, p = q = x = Q, ^ = —3) 
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Figure 22: Soliton on six quadrants riirij < 0 (oi = 1, 02 = 1.5, 03 = 2, 
k = 0.3, p = q = x = ^ = Qi) 
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7 Discrete linear wave equation 

It is interesting to compare the theory of the homogeneous dKdV equation 
we discussed above with the discrete linear wave equation 

V 12 -V 1 -V 2 + V = 0 (19) 


on quad-graphs. 

Obviously, equation (^) always can be solved with respect to any vertex 
and it does not contain singularities in contrast to the dKdV case. Therefore 
the splitting of self-intersecting strips is not valid in this case. 

The most serious distinction is that since the discrete wave equation does 
not contain the parameters the corresponding matrices (§) 



contain no spectral parameter and give the fake zero curvature represen¬ 
tation. In particular, the refactorization scheme from the Section ^ does 
not work, and the Theorem ^ is not valid, as the following simple exam¬ 
ples demonstrate. Of course, for the wave equation it is sufficient to analyze 
the propagation of delta-like initial data (black circle denotes 1, white one 
denotes —1, no circle denotes 0). 



a b c 


Figure 23: Interaction of the linear wave with the weak defect: splitting, 
reflection, bending 

Nevertheless we can claim that a weaker version of the theorem is valid 
also for the discrete wave equation. Let us consider the Cauchy problem for 
this equation on a regular square lattice with a weak defect D with initial 
data on a staircase (or coordinate half-axes) outside the defect. Again as in 
dKdV case according to our Theorem || the solution exists and unique and 
one can compare this solution with the one on the regular lattice without 
defects. 
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Let us consider the cross C{D) corresponding to the defect D as the 
union of all the characterstics passing through D. Besides the defect itself 
it contains also all the vertical and horizontal strips passing through the 
rectangle. 

Theorem 5. For the discrete wave equation a weak defect D does not affect 
the solution outside the cross C{D). 

The proof easily follows from the fact that for the discrete wave equation 
the difference v — vi is preserved when we are moving along the corresponding 
charactersitic. Actually the affected area is only a half of this cross (“thick 
hook”; see the picture below). For a general equation one should expect the 
whole quadrant bounded by this hook to be affected. 



Figure 24: The area affected by the defect D for the discrete linear wave 
equation 


8 Concluding remarks 

We have shown that for the integrable equations on the quad-graphs it is 
possible to formulate a criterion for the correctness of the Cauchy problems. 
The 3D consistency condition plays a crucial role in our analysis. 

The question is whether such a criterion exists for a general eqaution on 
quad-graphs. One can easily produce a lot of examples of graphs and IVP 
which are not sensitive to a particular form of the equation, but to describe 
all of them in some geometric way seems to be a hard problem. 

We have not discussed also global geometrical aspects (e.g. the quad- 
graphs on the surfaces of genus g). For some interesting results in this di¬ 
rection in the linear case we refer to Novikov and Dynnikov papers [13 - 
16]. 
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